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Abstract 

We analyze coherent two-color photoassociation of a Bose-Einstein conden- 
sate, focusing on stimulated Raman adiabatic passage (STIRAP) in free- 
bound-bound transitions from atoms to molecules. This problem is of par- 
ticular interest since STIRAP is predicted to be necessarily absent in the 
nondegenerate case [Javanainen and Mackie, Phys. Rev. A 58, R789 (1998)]. 
However, Bose-stimulation enhances the free-bound dipole matrix element for 
an atomic condensate, and photoassociative STIRAP turns out to be a viable 
mechanism for converting an atomic condensate to a molecular condensate 

with near-unit efficiency 
32.80.-t,32.80.Wr,03.75.Fi 
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As the field of Bose-Einstein condensation (BEC) in dilute gases continues to develop, 
candidate systems for BEC are being discovered on a regular basis. Besides the original 
alkalis there is the recent observation of BEC in hydrogen |J. Meanwhile, a spin- 
polarized helium condensate is not out of the question [|J, and recent theoretical work has 
shown that coherent photoassociation might be used to produce a degenerate molecular 
gas (MBEC) from an already-Bose-condensed sample of atoms ||,||. Here we consider the 
formation of MBEC using coherent photoassociation. 

Photoassociation (PA) is the process in which a pair of atoms interacts with a photon, 
thereby making a transition from the two-atom continuum to a bound state of the molecule. 
Quantization of the molecular dissociation continuum allows one to describe such free-bound 
transitions using the standard techniques of few- level quantum optics ]7|,|8[]. Beyond providing 
a few-level framework, such a quasicontinuum approach has also led to a matter-quantized 
formulation of photoassociation analogous to the analysis of second-harmonic generation 
of light HI!. In this context of nonlinear matter optics ||10|| , atom-molecule conversion 
displays coherent BEC-MBEC oscillations, adiabatic following, and non-classical collapse 
and revivals PJTT|. 

Nonetheless, one-color free-bound photoassociation generally occurs to an excited elec- 
tronic state of the molecule, and the subsequent irreversible losses, whether due to pho- 
todissociation (PD) or spontaneous decay, tend to negate the benefits of the coherence. 
In this Letter we develop coherent two-color free-bound-bound photoassociation, where the 
primary photoassociated molecules are transferred with another laser field to a stable molec- 
ular state. The goal is to state-selectively convert a BEC into an MBEC, while averting 
irreversible losses from electronically excited molecules. We consider pulsed free-bound and 
bound-bound couplings that occur in the counter-intuitive order ||12|| , and so correspond 
to stimulated Raman adiabatic passage (STIRAP) from atoms to molecules. Beyond its 
distinctly nonlinear character, this problem is of particular interest since we have previously 
argued for the absence of free-bound-bound STIRAP J7|. However, the present work illus- 
trates that such transitions are possible in the case of a condensate, since the free-bound 
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dipole matrix element is Bose-enhanced due to the fact that many atoms are in the same 
quantum state. A similar statistical effect may also increase hot pion production in nuclear 
physics |Uj , and speed up the decay rate of the inflation field in cosmology ]T1| . 

The development herein is outlined as follows. First, we review the many-body en- 
hancement for the PA dipole matrix element of a condensate, and discuss the implications 
for coherent free-bound-bound STIRAP. Next we take a semi-classical approach to identify 



the nonlinear counterpart of the dark state |15[ that contains no primary photoassociated 



molecules. Given that Bose-enhancement allows a counter-intuitive pulse scheme, the limit 
t — ► — oo finds the dark state with all atoms, while in the limit t — > +00 the dark state 
contains only molecules. It then appears possible that free-bound-bound STIRAP of an 
atomic condensate will produce a similarly degenerate gas of molecules. The remaining 
work establishes conditions that allow such adiabatic following to occur. 

Turning to the situation of Fig. [I], we assume that N identical atoms of mass m have 
condensed into the same one-particle state, say, a plane wave state with wave vector k = 0. 
Photoassociation removes two atoms from this state |1) and creates an excited molecule in 
state 1 2), with the translational energy due to the recoil momentum of the free-bound photon 
equal to % c^/Am. Including a second laser frequency, bound-bound transitions remove 
excited molecules from state |2) and create stable molecules in state |3). The translational 
energy per stable molecule arises from the combined recoil momenta of two photons, and 
equals fr 2 (qi — q2) 2 /4m. In second- quantized notation, we denote the Boson annihilation 
operators for atoms, primarily photoassociated molecules, and stable molecules, respectively, 
by a, b, and g. 

The laser-matter interactions that drive the atom-molecule and molecule-molecule tran- 
sitions are written in terms of their respective Rabi frequencies, k = d\E]_/2h and Q = 
d2E 2 /2h . Here the amplitude of the electric field driving a given transition is , and di 
is the corresponding dipole matrix element (i = 1,2). Lastly, we define the two-photon 
and intermediate detunings, A and S, both with due inclusion of photon recoil energies. 
Analogously to Ref. ||, the Hamiltonian for the system is 



— = fAa+a + Stfb - \k (aatf + aWfe) - Q (V + tfg) . (1) 

The Bose-enhancement of the free-bound dipole matrix element demonstrated as follows. 
First we consider the Heisenberg equations of motion, which determine the time evolution 
of the system according to 

a = —i (^2^ a ~ Ka ^^j 5 (2a) 
b = —i (sb — -kclcl — Qg \ , (2b) 
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g = ittb. (2c) 
Now, since the number of particles is conserved, 

0+0 + 2(6+6 + ^) =JV, (3) 



it is clear that a,b,g ~ yN. Hence, we define scaled boson operators of order unity as 
x — > x' = x/\/N, with x = a,b,g. Dropping the primes, the conserved quantity ([3]) is 
normalized to unity, and the equations of motion are given by 

a = -i (-Aa - x^bj , (4a) 
b = -i (db - ^xaa - Qg^j , (4b) 
g = iVtb . (4c) 

From Eqs. (|J), the many-body Bose-enhancement of the free-bound dipole matrix element 
is evident in the scaled Rabi frequency x — K - 

As it happens, the present STIRAP analysis depends crucially on the fact that the bare 
free-bound coupling k is scaled by the factor a/ZV, while that for the bound-bound transition, 
Q, is unchanged. To see why, we recall photoassociation in terms of our quasicontinuum 
model which represents the dissociation continuum of the molecule as an infinite num- 
ber of discrete energy levels separated equally by he. The correct free-atom results are 
recovered by taking the limit e — > 0. But, the bare free-bound coupling is given in terms 
of the frequency spacing as k ~ y/e, or in terms of the corresponding quantization volume 



as K ~ 1/y/V. In the free-atom limit (e — > or V — ► oo) we then find that k — ► 0. While 
this observation does not condemn photoassociation of a nondegenerate gas in the thermo- 
dynamic limit J7j^|, V — >• oo and A" — > oo with p = A^/V constant, it does imply that 
K(t) <C O(t) for all t. Hence, a counter-intuitive reversal of the coupling strengths of the 
pulses ||12|| , Q(t) ^> K,(t) for t — > — oo and /c(i) ^> fi(t) for t — > +00, cannot be achieved. It 
is therefore the absence of STIRAP that has been predicted for nondegenerate free-bound- 
bound transitions f7|. However, as shown above, in a condensate the bare Rabi frequency 
of a nondegenerate gas k is scaled by the Bose enhancement factor \^N, which leads to a 



finite value even in the thermodynamic limit ||; % ~ yN/V ~ yfp. This observation will 
open the door to using STIRAP as a means to create a stable molecular condensate. 

In order to facilitate an analytical solution, we define the "Kamiltonian" for the system 
by adding a multiple of the conserved particle number to the Hamiltonian 

K = H -hn\a j a + 2(tfb + g j g)} , (5) 



where the real constant h/j, is identified as the chemical potential per atom. The Heisenberg 
equations of motion for the unit-scaled operators become 
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— -A ) a + x a ^b 



b = i 



(2/i - 5) b + ^xaa + &g 



(6a) 
(6b) 



g = i (2fig + Sib) . (6c) 

From this point onward we also resort to the semi-classical approach analogous to the Gross- 
Pitaevskii approximation used to describe an alkali condensate. Accordingly, the quantities 
a, b and g are now c-numbers, rather than operators. 

We are looking for adiabatic solutions to the time-evolution equations (0) for transient 
couplings x(t) an d Q(t). Denoting the characteristic Rabi frequency scale for the light pulses 
by R and the characteristic pulse width by T, the adiabatic approximation (x ~ 0, x — a, b, g) 
should be valid when the evolution time scale of the system is short compared to the time 
scale of the pulses; for instance when RT 3> 1. We return to the adiabatic condition in a 



moment, and for the time being simply assume time scales for the problem such that x ~ 
is valid. 

Besides the steady state, we now specify exact two-photon resonance (A = 0). To 
economize the ensuing expressions, we also choose the Bose-enhanced Rabi coupling \ as 
the scale for our frequencies by writing Q = fix, 5 = 5x- Now, we have already assumed the 
Rabi frequencies to be real. Without sacrificing generality, this will allow us to consider the 
amplitudes a, b, and g as strictly real. Additionally, for counter-intuitive pulses the limits 
t — > (— oo, +oo) correspond to the limits Q — > (oo, 0). Thus we neglect any solution which 
is not similarly real for all values of fl from to oo. Discarding also solutions differing only 
by redundant signs, there simply remains 

fio = 0, (7a) 
a = ^Cl ^^2 + 0? - , (7b) 
bo = 0, (7c) 



«7o = -- (V2 + ^ 2 -fi); (7d) 



n 2 + 5 (5 ± Vfi 2 + 6 2 
/i ± = ± -^=— L x, 



2V^ 2 + 5 2 



± 2 



\ 



1 ± 



X 2 + 5 



5 
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(8a) 
(8b) 

(8c) 



9± = —x* 



1 ± 



2\ VCl 2 + 5 2 



n 2 + 5 -5 ± Vn 2 + 5 2 



nVn 2 + 5 2 



(8d) 



It is easy to see that the above results are the nonlinear counterparts of those obtained 
from the standard STIRAP analysis []12[ . The solution (|?p is the dark state. For a counter- 
intuitive pulse sequence, Cl ^> 1 for t — > — oo and Cl <ti 1 for t —>■ +oo, the dark state (0) 
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initially consists of atoms (ciq — > 1, <7o — > 0), while the final dark state is all molecules (do — > 
0, po ~~ > l/v2)- At no intervening time is there any population in the intermediate molecular 
state (&o = 0). If the laser pulses allow for adiabatic evolution, an atomic condensate is 
converted to a molecular condensate without any loss from the intermediate state. 

We have checked these results with an exact numerical solution to the equations of mo- 
tion (ID, examining the intuitive adiabatic condition RT > 1 as well. In Fig. |2|, Gaussian 
pulses of the form X (t) = Rexp [-(* - L>i) 2 /T 2 ] and Q(t) = Rexp [-(* - D 2 ) 2 /T 2 } illus- 
trate that an insufficient pulse area (RT) results in a significantly populated intermediate 
state. On the other hand, a pulse area much greater than unity readily decouples the nona- 
diabatic states @, allowing the system to adiabatically follow the dark state as it moves 
from the initial BEC to MBEC. In particular, Fig. 0(b) gives the probability of creating 
excited molecules as \b\ 2 ~ 10~ 7 , so that irreversible losses from either photodissociation 
or spontaneous decay should be negligible. Coherent free-bound-bound STIRAP is thereby 
confirmed. 

To improve upon our discussion of the adiabatic approximation, we focus specifically 
on the effect of pulse overlap. We apply the textbook criterion |16| for adiabaticity to the 
eigenvector if) = (a,b,g) T . From Eqs. (|7|) and (|8]), the coupling between the nonadiabatic 
states, i[)±, and the rate of change of the adiabatic state, ipo, must therefore be much less than 
the spacing between the respective chemical potentials, ip±ipo <C \fio — f/,±\. Furthermore, 
we restrict our analysis to zero intermediate detuning (5 — 0), and introduce the arbitrary 
pulse shapes x(t) = Rfi( T ) an d £l(t) = Rf 2 (r), where r = t/T is a dimensionless time. The 
adiabatic condition is then F <C RT, where the (dimensionless) nonlinear adiabatic factor, 
F, is given as 

F \h(r)lh(r)-h(r)lf 2 (T)\ 

2pft(T)+ft(T) 

x / 2 (r)//i(r) - y/2 + fi(T)/f?(T) | . (9) 

We find a dependence on pulse overlap that is qualitatively similar to the results in 
Ref. ]12| for ordinary STIRAP. In particular, if the two pulses vanish concurrently (such that 



their ratio is finite), the fraction in Eq. (^) diverges, and the adiabatic condition, F <C RT, 
is violated. Additionally, an increase in pulse area RT may provide adiabaticity despite 
a poor pulse overlap. Using Gaussian pulse shapes, these observations may be quantified 
as follows. In the vicinity of pulse overlap, at dimensionless times r ~ (D\ + D-ijjlT, we 
determine the maximum value of F(t), -F max , as a function of the delay between the pulses, 
D = D\ — D2- Numerically, breakdown of adiabaticity occurs at a specific pulse separation 
when the fractional efficiency, 2\g\ 2 , is no longer of order unity; hence, the value of F max at 
this point defines the adiabatic condition in terms of D. The results for RT = (10 2 , 10 3 , 10 4 ) 
are shown in Fig. |3|. We find that the adiabatic approximation is valid for pulse overlaps 
satisfying F max (.D) ~ 0.25 RT. Incidentally, we find exactly the same result for peaked- 
exponential pulse shapes, fi^ij) = exp(— |t — D 12 \/T). 

It remains to discuss a few items that we have so far ignored. The first regards the 
explicit role of photodissociation (PD). Consistent with the quasicontinuum model [[7|.|8|.pl] . 
we consider PD as irreversible decay from the intermediate state to the atomic continuum. 
Hence, if there is no intermediate-state population, there is no photodissociation. Second, 
our results are of course valid regardless of trapping of the atoms and molecules, provided 
that the time scale for coherent free-bound-bound STIRAP is shorter than the time scale 
for the motion of the atoms and/or molecules in the trap. Third, if laser intensities permit 
STIRAP during a time much shorter than the time scales for collisions between atoms and 
molecules, collisions are negligible as well. Fourth, we have not mentioned the divergence of 



the adiabatic factor in the wings of the pulses [12|. This is an artifact due to the mathemat 



ical shape of the pulse, and can be avoided if, in the region of interest, the shapes satisfy 
d(\ogQ) /dr = 0. In fact, this is the case for the peaked exponentials mentioned above. 

In conclusion, we hold the line on the absence of STIRAP in non-degenerate free-bound- 
bound transitions, while at the same time proposing the counter-intuitive pulse scheme as a 
possible mechanism for creating a molecular condensate from an initial BEC. This dichotomy 
arises because, for the case of a condensate, all N atoms are in the same quantum state, and 
the subsequent Bose-enhancement of the free-bound dipole matrix element enables a counter- 
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intuitive reversal of the Rabi frequencies associated with the light pulses. Our numerical 
trials have confirmed that STIRAP should take place in two-color photoassociation, and 
validated a simple quantitative criterion for adiabatic atom-molecule conversion. 

This work was supported in part by the NSF, Grant No. PHY-9801888, and the Research 
Experience for Undergraduates Program. Additional support provided by NASA, Grant No. 
NAG8-1428 and the Connecticut Space Grant College Consortium. 
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FIGURES 




FIG. 1. Three-level illustration of coherent free-bound-bound photoassociation, where N atoms 
have assumedly Bose-condensed into state |1). The free-bound and bound-bound Rabi frequencies 
are k and Q, respectively. Similarly, the two-photon and intermediate detunings are A and S. 
Irreversible losses are due to either photodissociation, at the rate T, or spontaneous decay, at the 
rate r$- The difference from the familiar three- level scheme is that the free-bound interaction 
involves three particles, and is therefore nonlinear. 
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FIG. 2. Transient free-bound-bound photoassociation of a BEC for a counterintuitive pulse 
order. The pulses are of equal height, xo = = R, so that R = 1 sets the unit of frequency. 
The intermediate detuning is 5 = 1, and the pulse delays are Z?i = 4.5 T and Z?2 = 2.5 T. Also, 
recall that ./V atoms gives a maximum of N/2 molecules, so that \g\ 2 = 1/2 represents complete 
conversion of the initial condensate, (a) T = 10 results in a pulse area insufficient to provide 
adiabaticity. (b) As expected, at T = 10 4 the initial BEC is converted entirely into MBEC. 



12 




2 4 6 8 10 

D/T 



FIG. 3. Efficiency of free-bound-bound STIRAP and the nonlinear adiabatic factor F plotted as 
a function of the pulse separation D = D\ — D2. In the region of pulse overlap, r « (D\ + Z?2)/2T, 
we determine the maximum value of F(t), F max , and compare it to the numerical results for the 
fractional efficiency, 2\g\ 2 . For R = 1 and T = (10 2 , 10 3 , 10 4 ), it is clear that F max /(0. 25 RT) = 1 
exactly marks the breakdown of adiabaticity. 
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